Evolution from BCS superconductivity to Bose-Einstein condensation: Current 
correlation function in the broken-symmetry phase 
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We consider the current correlation function for a three-dimensional system of fermions embedded in 
a homogeneous background and mutually interacting via an attractive short-range potential, below 
the (superconducting) critical temperature. Diagrammatic contributions in the broken-symmetry 
phase are identified, that yield for the (wave-vector and frequency dependent) current correlation 
function the fermionic BCS approximation in the weak-coupling limit and the bosonic Bogoliubov 
approximation in the strong-coupling limit (whereby composite bosons form as bound-fermion pairs) . 
The temperature dependence of the superfluid density (from the BCS exponential behavior at weak 
coupling to a power-law behavior at strong coupling) and the form of the Pippard-like kernel at 
zero temperature are explicitly obtained from weak to strong coupling. Quite generally, it is shown 
that the Pippard-like kernel is the sum of a local (London-like) term and of a nonlocal compo- 
nent, the local term being dominant in the strong-coupling limit and the nonlocal component in 
the BCS (weak-coupling) limit. It is also shown that the range of the nonlocal component is de- 
termined by the coherence length measuring the spatial correlations of the amplitude of the order 
parameter, namely, the correlations among different Cooper pairs (or composite bosons), rather 
than between the fermionic partners of a given pair. In addition, a prescription is provided for 
mapping the fermionic onto the bosonic diagrammatic theories in the broken-symmetry phase, thus 
complementing what already done in the normal phase. 



PACS numbers: 74.20.-z,74.20.Fg,74.25.Nf 
I. INTRODUCTION 

Already from London and Pippard phenomenological 
approaches to superconductivity, the current response to 
an external magnetic field played a key role in account- 
ing for the Meissner effect. The Meissner effect was, in 
fact, later demonstrated explicitly within the microscopic 
BCS approach, by examining the behavior of the (static) 
current correlation function 1 ' 2 . 

The BCS approach as it stands applies to weak- 
coupling superconductors, with largely overlapping 
Cooper pairs. Recently, however, it has become of inter- 
est to study the evolution of a superconducting system 
from weak to strong coupling, with the strong-coupling 
limit corresponding to a system of nonoverlapping (com- 
posite) bosons 3 " 9 . The physical motivation to these 
works stems essentially from cuprate superconductors, 
whereby the small value of the (superconducting) coher- 
ence length as well as the presence of a pseudogap above 
the critical temperature in the underdoped region ,n 
have suggested the possible relevance of a crossover sce- 
nario from a weak-coupling limit for overdoped samples 
(with Cooper pairs forming and condensing at the critical 
temperature within a BCS description) toward a strong- 
coupling limit for underdoped samples (with preformed 
(composite) bosons existing above the superconducting 
critical temperature and Bose-Einstein (BE) condensing 



below it). For cuprate superconductors, it seems actu- 
ally likely that the pairing corresponds to an intermediate 
regime between overlapping Cooper pairs and nonover- 
lapping composite bosons. More generally, the possible 
occurrence of a crossover from BCS to BE has been em- 
phasized for other classes of "unconventional" supercon- 
ductors as well 12 . 

The first discussion of the crossover from BCS to 
BE can be found in the pioneering paper by Ea- 
gles dealing with superconductivity in low-carrier doped 
superconductors 13 . A systematic approach to this prob- 
lem was then provided by Leggett 14 , who showed that 
for a fermionic system with an attractive interaction at 
zero temperature a smooth crossover from a BCS ground 
state of overlapping Cooper pairs to a condensate of com- 
posite bosons occurs as the strength of the attraction 
increases. A few years later the problem was reconsid- 
ered by Nozieres and Schmitt-Rink 15 , who extended the 
approach to finite temperatures by using diagrammatic 
methods, always finding a smooth crossover between the 
two limits. 

Several quantities have so far been calculated within 
the BCS-BE crossover approach, including the criti- 
cal temperature and chemical potential 16 , the coherence 
length(s) at zero temperature 6 , the spectral function 
above 17 and below 18 the critical temperature, and the 
effect of an external magnetic field on the pseudogap tem- 
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perature 19 . No systematic study yet exists, however, for 
the current correlation function and derived quantities, 
such as the (temperature dependence of the) superfiuid 
density and the Pippard kernel. Purpose of this paper is 
to provide a detailed calculation of these quantities, by 
relying on controlled approximations in the two (weak- 
and strong-coupling) limits, which allow us to examine 
specifically how the response of the original Fermi sys- 
tem can be interpreted in terms of the response of an 
effective Bose system in the strong-coupling limit. 

We shall rely on a diagrammatic approach that selects 
the relevant contributions to the (wave-vector and fre- 
quency dependent) current correlation function in the 
weak- and strong-coupling limits separately, and then 
sum these contributions to obtain the current correlation 
function for all couplings. The contributions selected in 
the strong-coupling limit will appear as fluctuation con- 
tributions when extrapolated to the weak-coupling limit. 
In this paper, we will consider a Fermi system with an at- 
tractive (point-contact) interparticle interaction embed- 
ded in a three-dimensional continuum, with an isotropic 
(s-wave) gap function. The ensuing diagrammatic the- 
ory in the broken- symmetry phase will be thus simpli- 
fied considerably, yet reproducing the desired features in 
both limits, in a similar fashion to what was done in the 
normal phase above the critical temperature 9,20 . No lat- 
tice effects nor the occurrence of an angular dependent 
(d-w&ve) gap function will accordingly be taken into ac- 
count, deferring their study to future work. The present 
approach should be regarded as a necessary preliminary 
step toward the development of a theory of the current 
response function over the whole crossover range for more 
realistic models. 

As far as the current correlation function of inter- 
est is concerned, the two (weak- and strong-coupling) 
limits will be described by the BCS approximation 
(for fermions) and the Bogoliubov approximation (for 
bosons), respectively. Again, these represent the sim- 
plest approximations which can be conceived to de- 
scribe fermionic and bosonic systems, respectively, in 
the broken-symmetry (superfiuid) phase, and share anal- 
ogous self-consistent constructions to generate the rel- 
evant excited states within a temperature-dependent 
mean-field approximation. Albeit the Bogoliubov ap- 
proximation for bosons is known to suffer from several 
shortcomings 21,22 (which could be overcome, at least in 
principle, by considering more involved approximations 
for condensed bosons 23 ), we shall anyhow restrict our- 
selves to this approximation for the bosonic limit, since 
we are ultimately interested in obtaining a reasonable 
strong-coupling limit while following the evolution of 
the fermionic current correlation function from weak to 
strong coupling. In this respect, we should mention that 
the Bogoliubov approximation has been recently adopted 
to get the response function near zero temperature for a 
Bose gas in the context of the BE condensation of ultra- 
cold atomic dilute gases 24 . 

While the BCS approximation to the current corre- 



lation function is well established in weak coupling 1,2 , 
recovering the Bogoliubov results for (composite) bosons 
from an originally fermionic theory is altogether nontriv- 
ial and will therefore be considered at some length in 
this paper. Specifically, we shall identify the fermionic 
diagrams (over and above the standard BCS "bubble" 
contribution) which reproduce both the longitudinal and 
transverse contributions to the (wave-vector and fre- 
quency dependent) Bogoliubov current correlation func- 
tion, namely, the "ladder" diagrams known for making 
the BCS result conserving 1 and the Aslamazov-Larkin 
(AL) type diagram(s) 25 here extended to the broken- 
symmetry phase, in the order. Since these diagrams can 
be evaluated analytically in the strong-coupling limit, a 
detailed account of this derivation will be provided in 
what follows. In this context, the prescription for map- 
ping the fermionic onto the bosonic diagrammatic theo- 
ries in the broken-symmetry phase will also be provided, 
in an analogous fashion to what already achieved in the 
normal phase 9 . 

Besides providing a formal theory that bridges weak- 
and strong-coupling approximations to the current cor- 
relation function, the main physical results of this paper 
concern, in particular, the temperature dependence of 
the superfiuid density (from the BCS exponential behav- 
ior at weak coupling to a power-law behavior at strong 
coupling) and the form of the Pippard-like kernel at zero 
temperature 26 . The contributions of the BCS and AL di- 
agrams to the temperature dependence of the superfiuid 
density will be numerically calculated for various cou- 
plings from weak to strong, and their relative importance 
assessed. [These calculations require only the static con- 
tribution to the current correlation function, so that the 
analytic continuation from Matsubara to real frequencies 
is not required.] It will be also shown that the Pippard- 
like kernel is the sum of a local (London-like) term and of 
a nonlocal component, where only the local term survives 
in the strong-coupling limit, thus correctly reproducing 
the result for non-interacting (composite) bosons, while 
only the nonlocal component survives in the BCS (weak- 
coupling) limit, as expected. It will be further shown that 
the characteristic length that controls the range of the 
nonlocal component coincides (for a clean system) with 
the coherence length £ p hase introduced in Ref. 6 as a mea- 
sure of the spatial correlations of the amplitude of the or- 
der parameter, over the whole crossover range from weak 
to strong coupling. The length £ p hase relates the corre- 
lations among different Cooper pairs or between com- 
posite bosons, while correlations between fermionic part- 
ners in a pair is associated with the different lenght scale 
Cpair- 6 The quantity £ p hasc, in turn, coincides with the 
Pippard coherence length £ (and with £ pa ir) in the weak- 
coupling limit, decreases as the coupling is increased from 
the weak- toward the intermediate-coupling regime, and 
eventually starts again to increase upon approaching the 
strong-coupling (bosonic) limit where £ p hase ^ £pair 6 (see 
also Ref. 27 for an extension to the lattice case). 

The plan of the paper is as follows. In Section II we 



2 



consider the current correlation function within the (con- 
serving) BCS approximation (i.e., bubble plus ladder di- 
agrams) and show that, in the strong-coupling limit, it 
maps onto the longitudinal part of the current correla- 
tion function for (composite) bosons within the Bogoli- 
ubov approximation. In Section III we identify, in addi- 
tion, the contribution to the fermionic current correlation 
function that maps onto the transverse part of the cur- 
rent correlation function for (composite) bosons within 
the Bogoliubov approximation. In Section IV numeri- 
cal results for the superfluid density vs temperature as 
well as for the Pippard-like kernel at zero temperature 
are presented, over the whole range from weak to strong 
coupling. Section V gives our conclusions. The Appen- 
dices discuss more technical material. Specifically, in Ap- 
pendix A the connection between the particle-particle 
ladder in the broken-symmetry phase and the Bogoli- 
ubov propagators is obtained; in Appendix B the map- 
ping between the fermionic and bosonic diagrammatic 
structures in the broken-symmetry phase is established; 
in Appendix C the form of the wave- vector and frequency 
dependent current correlation function for bosons within 
the Bogoliubov approximation is discussed, with empha- 
sis on the ensuing form of the Pippard-like kernel; and in 
Appendix D the BCS bubble is examined in a novel way 
as to extract its kernel for all couplings and not only in 
the weak-coupling limit. 



II. CURRENT CORRELATION FUNCTION 
WITHIN THE BCS APPROXIMATION IN THE 
STRONG-COUPLING LIMIT 

In this Section, we reconsider the well-known BCS 
approximation for the current correlation function 1 , in 
the form which fulfills conservation laws 28 . We will 
specifically show that in the strong- coupling limit of the 
fermionic attraction (and at temperatures small com- 
pared with the dissociation energy of composite bosons) 
this approximation maps onto the longitudinal part of 
the current correlation function for bosons within the 
Bogoliubov approximation. This result, which is per se 
nontrivial, requires us to search for additional contribu- 
tions to the current correlation function which, in the 
strong-coupling limit, reproduce the transverse part of 
the current correlation function for bosons within the Bo- 
goliubov approximation, as discussed in the next Section. 

Quite generally, the current correlation function at fi- 
nite temperature in the imaginary time (Matsubara) rep- 
resentation reads 2 : 



X 7l7 '(rT,rV) = -<T r [j 7 (rr)j 7 ,(rV)]) 



(1) 



where < r, r' < (3 ((3 — l/(fcsT) being the inverse 
temperature and fc^ the Boltzmann's constant), T T is 
the imaginary-time ordering operator, (• • •) is a thermal 
average, and 7, 7' = (x, y, z) label Cartesian components. 
For our purposes, it is convenient to write the current 



operator in (1) in terms of the Nambu representation of 
the field operators 
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as follows: 



j« = 2L E (^( r ) V ^( r ) - (V*J(r))*/(r)) (2) 

m being the fermionic mass. 

The current correlation function (1) can be expressed 
in terms of the single-particle Green's function 

a(l,2) = -<T T [*(l)*t( 2 )]) (3) 
and of the two-particle Green's function 

2 (1, 2; 1', 2') = (T T [*(l)*(2)*t(2')* t (l')]> (4) 

(with the short-hand notation 1 = (ri,n,^i) and so on, 
in terms of the Nambu spinor components) as follows. 
Let 

L(l, 2; 1', 2') = G 2 (l, 2; 1', 2') - 0(1, l')0(2, 2') (5) 

be the two-particle correlation function, which satisfies 
the Bethe-Salpeter equation 

L(l, 2; 1', 2') = -0(1, 2')0(2, 1') + J d3456 0(1, 3) 

x 0(6,l')5(3,5;6,4)L(4,2;5,2') (6) 

where 



B(1 2-1' 2') * S(M/) 
4 " ' J " 50(2', 2) 



(7) 



is an effective two-particle interaction. Equation (6) 
can be formally solved in terms of the many-particle T- 
matrix, defined as the solution to the equation 29 

T(l, 2; 1', 2') = 3(1, 2; l', 2') + J rf3456 3(1, 4; l', 3) 

x 0(3,6)0(5,4)T(6,2;5,2') (8) 

by writing 

-£(l,2;l',2') =0(l,2')0(2,l') + j d3456 0(l,3) 

x 0(6, l')T(3, 5; 6, 4)0(4, 2')0(2, 5). (9) 

In terms of the above quantities, one obtains the following 
expression for the current correlation function (1): 

X 7 ,y(rr,rV) = (v 7 - V 7 ) (v' y - V?,) £ 

v ; 1,1' 

x L(r"r^,r"Vr;rT + ^,rV'+^)|r»=r,r'»=r'. 

(10) 
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In Fourier space, the above relation reads (in three di- 
mensions) : 



X 7> y (q, 00 - {2mi)2 g J ( J )3 | (2 P )3 ^ E • 

x ^£e^(2p 7 + g 7 )(2p;, + g 7 ,) 

n' 

x T«,(pw„,pW;q^) 



(11) 



where q, p, and p' are wave vectors, £2„ = 2-jtv//3 (u 
integer) is a bosonic Matsubara frequency, while io n = 
(2n + l)7r//3 and w n = (2n' + l)7r//5 (n, n' integers) are 
fermionic Matsubara frequencies (77 being a positive in- 
finitesimal) . The conventions for the (incoming and out- 
going) four- vectors (with q = (q, fi„), p = (p, w„), and 
p' = (p',c<j n ')) and for the Nambu indices are shown in 
Fig. 1, both for L and the many-particle T-matrix. 30 



p+q 




P l\ 



FIG. 1. Graphical representation of the two-particle cor- 
relation function (L) in terms of the many-particle T-matrix 
(T) (four- momenta and Nambu indices are indicated) . 

The above equations hold quite generally, regardless 
of the specific approximation for the kernel S defined in 
Eq. (7). The BCS approximation is made manifest at 
this point by setting: 



S B cs(l,2;l',2') 



5T, 



BCS 



(1,1') 



(12) 



SQbcs{2',2) 

= - T le 2 , S ( x i - x 2 ')V{x+ - x v ) 

x S(x v - x 2 )Tf il( , 2 (l - 5 eiil ,) 

where, as usual, only the off-diagonal terms of the BCS 
self-energy have been retained (r 3 being a Pauli matrix, 
and with the additional short-hand notation x\ = (iq, n) 
and so on). Note that the factor V(x 1 — xy) = S(t^ — 
Ty)V{rx — Ty) represents the (instantaneous) attractive 
effective fermionic interaction, as specified in Appendix 
A. 

The first term on the right-hand side of Eq. (9), when 
inserted in Eq. (10), yields the standard BCS "bubble" 
contribution to the current correlation function 1,2 , pro- 
vided the BCS single-particle Green's functions 

£(k) + iu> n 



£/n(k,w n ) 
£?2i(k,w„) 



A _ 



= -^22(-k, -U) n ) 

Gi 2 (k,u> n ) 



(13) 



are used. In the above expressions, £(k) = k 2 / (2m)— fi {ji 
being the chemical potential) and E(k) = ^/£(k) 2 + A 2 
for an isotropic (s-wave) gap function A. The BCS 
bubble contribution to the current correlation function 
taken alone does not fulfill the conservation laws (and, 
in particular, the longitudinal f-sum rule 1 ); yet it gives a 
reasonable account of the transverse current correlation 
function in the weak-coupling limit 1 ' 2 . 

Quite generally, the superfluid density p s can be ob- 
tained from the knowledge of the transverse component 
of the current correlation function as follows 31 : 



{n~ Ps) 



7,7' 



in 



(14) 



where n is the density and t is a unit vector transverse 
to q = q/|q|. In particular, inserting the BCS bub- 
ble contribution to the current correlation function in 
Eq. (14) produces a meaningful overall temperature de- 
pendence of the superfluid density in the weak-coupling 
limit, not only for an s-wave 2 but also for a <i-wave gap 
function 32,33 . In the strong-coupling limit (and for tem- 
peratures small compared with the dissociation energy 
of the composite bosons), on the other hand, the BCS 
bubble contribution to the current correlation function 
vanishes, thus requiring consideration of additional dia- 
grammatic contributions. 

Still within the BCS approximation, consideration of 
the second term on the right-hand side of Eq. (9) is re- 
quired to fulfill conservation laws, provided that BCS 
single-particle Green's functions (13) are consistently 
used 1 . In particular, it can be shown that the longi- 
tudinal f-sum rule 31 



X] *?7 X 7 ,7' (q -> °> = 0) qy 



n 
m 



(15) 



is fulfilled at any coupling in the static limit when this 
term is properly included. 

It may further be of interest to obtain the full wave- 
vector and frequency dependent current correlation func- 
tion when the kernel S is taken within the BCS approxi- 
mation. To this end, it is necessary to solve the integral 
equation (8) in Fourier space with the BCS form (12) 
of the kernel, obtained from the off-diagonal terms of 
the BCS self-energy. With this restriction, only four el- 
ements of the many particle T-matrix with £l ^ t' L and 
£r 7^ l' R survive (cf. Fig. 1) [we adopt the convention 
1 <-> (i = 1,1' = 2) and 2 <-> {t = 2,1' = 1) to label the 
nonvanishing matrix elements of the T-matrix] . In addi- 
tion, when the attractive effective fermionic interaction 
is taken of the form of a contact potential with strength 
v (cf. Appendix A), these four independent elements 
can be shown to satisfy the following algebraic equation: 



+ v 



Tu(q) T 12 (q) 
T 2 i{q) T 22 (q) 

-nn(«) 
n 2 i(<?) 



ni 2 (?) 
-iMg) 



1 
1 



Tu(q) 
T 2 i(g) 



Ti 2 (q) 
T 22 {q) 



(16) 
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In this expression 

dp 1 



nn (9) = J (||^ £ en (p + 9)Sn (-p) = n 22 (-«) 



(17) 



and 

ni 2 (g) = / ^3 ^£0i 2 (p + 9)012 (-p) = n 21 (g) 

(18) 

are the elementary rungs on which the many-particle T- 
matrix is built up within the BCS approximation. 

An equation similar to (16) holds for the particle- 
particle ladder in the broken-symmetry phase (cf. Ap- 
pendix A). Comparison of the two equations leads us to 
identify: 



Tn{q) = -T n (q), T 12 (q) = T 12 (q) 
T 21 (q)=T 21 (q), T 22 (q) = -T 22 (q). 



(19) 



In particular, the results obtained in Appendix A in the 
strong-coupling limit (whereby /3fi — ► — oo) can be ex- 
ploited to write: 



T u (q) = (8Tr/m 2 a F )gi(q) 
T 12 (q) = (8TT/m 2 a F )G 21 (q) 



(20) 



where ap is the fermionic scattering length (which 
is positive in the strong-coupling limit) and Q' and 
Q' 21 are the normal and anomalous noncondensate 
bosonic Green's functions, respectively, in the Bogoli- 
ubov approximation 2 . 

Within the above approximations, the many-particle 
T-matrix depends only on q and not on all the incoming 
and outgoing four- vectors (cf. Fig. 1). For this reason, 
the remainder (i.e., other than the bubble) contribution 
to the current correlation function within the BCS ap- 
proximation can be cast in the following form: 

X ™y(q) = I- y (q)T u (q)I ll (q)+I 7 (q)T 12 (q)I Y (-q) 

+ I 7 (-q)T 21 (q)Iy(q) + I 7 (-q)T 22 (q)I Y (-q) 

(21) 



where we have set 

T(\-l. f dp 1 

l[Q) ~ mj (27r)3/3 



£(2p 7 + q 1 )Qi\{p + 9)021 (p)- 



(22) 

It can be readily shown that in the strong-coupling limit 
I-y(q) becomes 



i-y{q) ^ -q-y 



majfA 
16tt 



(23) 



so that the expression (21) becomes eventually with the 
help of Eq. (20): 



(2m) 



2 [0'(9)+0'(-9)- 20 2 i(9)] 



(24) 



where the result (A37) (relating the general definition 
of the condensate density n — \a\ 2 to the off-diagonal 
single-particle BCS Green's function Q\ 2 in the strong- 
coupling limit) has been exploited. Note that, when the 
diagonal single-particle BCS Green's function On is used 
to evaluate the particle density n, a relation similar to 
(A37) holds for n/2 in the place of no, implying no = n/2 
at the level of the above approximations. 34 

The expression (24) is purely longitudinal. This ex- 
pression can be thus compared with the longitudinal part 
of the current correlation function for a system of bosons 
within the Bogoliubov approximation (cf. Appendix C). 
Identifying n/2 with the bosonic density n# and 2m with 
the bosonic mass tub (and neglecting the difference be- 
tween the bosonic density rig and the condensate den- 
sity no), one sees that the Bogoliubov result (C4) is four 
times smaller than the fermionic result (24) in the strong- 
coupling limit. This overall factor of four stems from hav- 
ing considered the current associated with the particle 
number and not with the charge density. [This is equiv- 
alent to setting formally e — 1 for the fermionic charge 
and es = 2 for the bosonic charge.] From Eq. (24) it 
can also be explicitly verified that the f-sum rule (15) is 
satisfied at this order of approximation. 

Since the current correlation function for bosons within 
the Bogoliubov approximation contains not only a longi- 
tudinal but also a transverse part (cf. Appendix C), in 
the strong-coupling limit the BCS approximation for the 
current correlation function clearly fails to recover the 
expression of the current correlation function within the 
Bogoliubov approximation for the composite bosons. It is 
therefore necessary to search for additional fermionic dia- 
grams for the current correlation function in the broken- 
symmetry phase, which recover the transverse Bogoli- 
ubov result in the strong-coupling limit. This task is 
achieved in the next Section. 



III. TRANSVERSE CURRENT CORRELATION 
FUNCTION IN THE BROKEN-SYMMETRY 
PHASE FROM THE ANALOGY WITH THE 
BOGOLIUBOV RESULT 

The simplest approximation for a Bose-condensed sys- 
tem is the so-called "Bogoliubov approximation", that 
relies on the system being "dilute" 35 ' 36 ' 2 . When one is 
interested in the crossover from a fermionic description 
of physical quantities in weak-coupling to a description 
in terms of (composite) bosons in strong-coupling, it is 
natural to require the (composite) bosons to be described 
by the Bogoliubov approximation, at least as a first sig- 
nificant step toward a more complete description of the 
strong-coupling limit. 
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Concernig specifically the current correlation function, 
the form (C4) for the other-than-longitudinal contribu- 
tion within the Bogoliubov approximation suggests that 
the fermionic diagrams to be considered should contain 
two particle-particle ladder propagators (of the type dis- 
cussed in Appendix A). In analogy with a similar connec- 
tion between fermionic and bosonic diagrams established 
in the normal phase 20 , the Aslamazov-Larkin 25 type dia- 
gram of Fig. 2(a) in the broken-symmetry phase appears 
to be a good candidate for our purposes. The correspond- 
ing contribution to the current correlation function reads: 

,(q\ = L_ MB_Iv f-^iv 

{W) (2m) 2 J (2tt) 3 ^ J (2tt) 3 fi 

x /t^3^EE E (2p 7 + 2g 7 + Q 7 ) 
X (2p' Y + 2qy + Q Y )T Vl {q)T vk (q + Q) 

x Qe L e L (p + q)Ge>[e L (p + q + Q)Ge[e L i~p) 

x Qi' R in ip' + q)Gt R t> R (p' + q + Q)Gt' R e» (-p') 

(25) 

where (p,p') and (q,Q) are fermionic and bosonic four- 
vectors, respectively, and the minus sign is due to the 
presence of an additional fermion loop with respect to 
the BCS fermionic bubble already considered in Section 
II. In the expression (25) we have taken into account 
the restrictions on the Nambu indices of the particle- 
particle ladder propagator, due to the regularization of 
the potential we have adopted. Note that all single- 
particle Green's functions in Eq. (25) are meant to be 
taken within the BCS approximation (13). 




P+q+' 




(b) 



p+q+Q 



FIG. 2. (a) AL type diagram for the current correlation 
function in the broken-symmetry phase; (b) Diagram topo- 
logically nonequivalent to (a) but with the same value for a 
contact potential. 

Although the expression (25) formally holds for all cou- 
pling, its calculation gets considerably simplified in the 
strong-coupling limit of interest, whereby the off-diagonal 



single-particle BCS Green's functions can be neglected 
with respect to the diagonal ones. In this case Eq. (25) 
reduces to: 



x# (Q)--f pjEE J'Hq.Q) r WL (q) 

xT iL e R (q + Q)J^(q,Q) (26) 
where we have introduced the vertex 
1 f dp 1 >— ^ 



J%Q) 



2mJ (2 , )3/3 ^( 2 ^ + 2 ^ + ^) 

X Gu{p + q)Gu(P + q + Q)Gu(-p) . (27) 



Since J 2 (g, Q) — J^{— q, —Q) owing to the symmetry re- 
lations in Eq. (13), it is sufficient to evaluate Jj(q,Q). 
One obtains: 



j'(q,Q) 



(2g 7 + Q 7 ) m 2 a F 



2m 



16tt 



(28) 



in agreement with the value obtained for the correspond- 
ing vertex of the AL diagram in the normal phase 20 . 
Making explicit the Nambu indices in Eq. (26) and re- 
calling Eqs. (19) and (20) for the connection between the 
elements of the particle-particle ladder and the bosonic 
propagators within the Bogoliubov approximation, we 
obtain eventually in the strong-coupling limit: 



1 1 



X7 < 7 '^ ~ (2m) 2 4 J (2tt) 3 



x (2q y + Qy) [g'(q)G'(q + Q) + G'(-q) 

x G'(-q - Q) - 2G' 21 (q)G' 21 (q + Q)\ ■ (29) 

This expression has to be compared with the current cor- 
relation function for bosons within the Bogoliubov ap- 
proximation, i.e., with the second term on the right-hand 
side of Eq. (C4). The comparison shows that expression 
(29) is only two times larger than its Bogoliubov coun- 
terpart, while one would expect it to be four times larger 
according to the discussion following Eq. (24) . The miss- 
ing factor of two is provided by an additional diagram 
[cf. Fig. 2(b)], which is topologically nonequivalent to 
the diagram of Fig. 2(a) but has the same value for a 
contact potential. 

In the weak-coupling limit, on the other hand, the 
AL type diagrams of Fig. 2 in the broken-symmetry 
phase represent fluctuation contributions to the ordinary 
BCS bubble approximation, and merge with the usual 
Aslamazov-Larkin contribution 25 when approaching the 
critical temperature from below. 

Having succeeded in obtaining the Bogoliubov result 
for the transverse part of the current correlation func- 
tion in the broken-symmetry phase as the strong-coupling 
limit of the fermionic diagrams of Fig. 2, we pass now to 
determine the evolution of the transverse part of the cur- 
rent correlation function from weak to strong coupling, 
by approximating the current correlation function in all 
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coupling regimes as the sum of the standard BCS bubble 
diagram and the diagrams of Fig. 2. In particular, we 
shall be interested in the static limit of these diagrams, 
from which the superfluid density p s can be extracted as 
q — > according to Eq. (14) and the Pippard kernel can 
be obtained by retaining finite values of q. The numer- 
ical results obtained, from weak to strong coupling, for 
the superfluid density as a function of temperature and 
for the Pippard kernel at zero temperature are discussed 
in the next Section. 



IV. NUMERICAL RESULTS FOR THE 
SUPERFLUID DENSITY AND THE 
PIPPARD-LIKE KERNEL FROM WEAK TO 
STRONG COUPLING 

In this Section, we evaluate numerically the AL type 
contribution (25) to the current correlation function for 
vanishing external frequency from weak to strong cou- 
pling, in addition to the standard BCS bubble contribu- 
tion. From this calculation we extract the temperature 
dependence of the superfluid density (when the external 
wave vector is also allowed to vanish) as well as the form 
of the Pippard-like kernel at zero temperature (by keep- 
ing the external wave vector finite). 

While the numerical calculation of the BCS bubble is 
rather straightforward, calculation of the expression (25) 
requires a non- negligible numerical effort. Specifically, 
evaluation of the particle-particle ladder T [cf. Eqs. (A9) 
and (A10)], where the frequency sum can be done ana- 
lytically, requires a two-dimensional numerical integra- 
tion; in addition, the four independent vertices with 
three single-particle Green's functions (of which the ex- 
pression (27) represents a degenerate case), where again 
the frequency sum can be done analytically, require a 
two-dimensional (three-dimensional) numerical integra- 
tion for vanishing (finite) external wave vector; finally, 
both the frequency sum and the wave vector integration 
over the bosonic variable q of the T functions have to be 
done numerically. Special care has to be exerted when 
performing this frequency sum, which requires explicit 
inclusion of about 300 to 500 (positive) Matsubara fre- 
quencies, as well as a careful treatment of the tail of 
the frequency sum which we approximate by an integral 
extending up to infinity. The resulting function (to be 
eventually integrated over the wave vector q) is, in fact, 
affected by spurious oscillations for large values of |q|, 
which result from incomplete cancellations of large num- 
bers when the tail of the above frequency sum is added to 
the contribution of the finite set of Matsubara frequen- 
cies. In addition, these oscillations are strongly ampli- 
fied by the factor q 4 originating from the measure of the 
(three-dimensional) integral and from the wave-vector 
dependence of the vertices with three Green's functions. 
As a consequence, we have pragmatically truncated the 
integral over |q| where these spurious oscillations start 



to appear. We have consistently verified that this trun- 
cation procedure yields correct results in the two cases 
where analytic controls can be made, namely, in the 
strong-coupling limit at low enough temperatures and in 
the weak-coupling limit close to the (BCS) critical tem- 
perature. 

A comment on the procedure we have adopted to gen- 
erate the values of the gap function A and chemical po- 
tential /i is in order at this point (these parameters enter 
the single-particle Green's functions (13) and, through 
them, all other relevant quantities). According to a stan- 
dard procedure when following the evolution from BCS 
superconductivity to BE condensation 13 ' 14 , A and /i are 
obtained by solving the coupled equations for the mean- 
field gap and particle number, in terms of the Green's 
functions (13). Strictly speaking, these equations are 
valid at low enough temperature (compared with the 
superconducting critical temperature), and corrections 
to these equations should accordingly be considered at 
higher temperatures. In this paper, we restrict to low 
enough temperatures, so that inclusion of these correc- 
tions is not expected to be important. Improved results 
at higher temperatures may then be obtained by sup- 
plying the present calculation for the current correlation 
function with improved values of A and fi, which take 
into account the effects of superconducting fluctuations 
below T c . 34 

A. Superfluid density 

Once the AL type and the fermionic bubble contribu- 
tions to the current correlation function are calculated 
as described above, the temperature dependence of the 
superfluid density p s can be obtained via the definition 
(14) by setting Q u = and letting q — > 0. 

The low-temperature behavior of the fermionic bubble 
leads to the well-known exponential suppression of the 
normal density n — p s when approaching zero tempera- 
ture; this exponential behavior is governed by the BCS 
gap in weak coupling and by the magnitude of the chem- 
ical potential in strong coupling. These energy scales, in 
turn, limit the respective temperature range for the va- 
lidity of the exponential behavior. The low-temperature 
behavior of the AL type diagram in the strong-coupling 
limit, on the other hand, yields a T 4 increase of the 
normal density for increasing temperature, in agreement 
with the Bogoliubov-Landau expression 2 (cf. Appendix 
C). In this limit, the T 4 behavior is valid for temperatures 
small compared with the bosonic chemical potential and 
is due to the linear dispersion of the Bogoliubov mode 
for small enough wave vectors. Since the poles of the 
particle-particle ladder T yield a linear dispersion even 
at intermediate and weak coupling 37 (giving origin to the 
so-called Bogoliubov- Anderson mode 38 ), one expects the 
T 4 behavior from the AL type contribution to survive at 
low enough temperature for all coupling ranges. This T A 
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behavior is thus expected to dominate over the BCS ex- 
ponential behavior at low enough temperatures. [Note, 
however, that the T behavior is bound to cross over into 
a different power-law behavior at a finite characteristic 
temperature, as discussed below.] 



1 . 4 




1 — 1 1 1 

-2-10 1 2 

(k F a F ) _1 



FIG. 3. Temperatures below which (dashed-dotted line) 
the AL type contribution to the normal density is dominant 
with respect to the BCS contribution and below which (dot- 
ted line) the T 4 behavior of the AL type contribution sets in, 
versus the coupling parameter (JifOf)" 1 . The superconduct- 
ing critical temperature T c (full line) is also shown for com- 
parison. [All temperatures are in units of the Bose-Einstein 
condensation temperature Tbe = 3.31n 2 J 3 /tub] 

An interesting question that can be answered with 
the present calculation is to determine the temperature 
range below which the AL type contribution dominates 
over the BCS contribution on the whole coupling range, 
from weak to strong. To this end, we have calculated 
the temperature at which these two contributions equal 
each other for several values of the coupling parameter 
(kpap) 1 (where kp is the Fermi wave vector and ap 
is the scattering length - cf. Appendix A). [In practice, 
the whole crossover region can be spanned by varying 
(kpap)^ 1 about from —1 to +1; deviations from weak- 
coupling approximations when (kpciF) 1 — 1 and from 
strong-coupling approximations when (kpap)^ 1 ^ 1 be- 
come, in fact, negligible for all practical purposes.] In 
Fig. 3 this temperature is plotted versus {kpap)" 1 , to- 
gether with the superconducting critical temperature T c 
obtained following Ref. 15 (cf. also Refs. 16 and 17). In 
the same plot, we have also reported the temperature 
at which the AL type contribution changes from the T 4 
behavior to a different power law, from which one notes 
that the range of validity of the "low-temperature" be- 
havior is significantly enhanced in the crossover region. 
[It is, however, clear from the previous discussion that 
the results for these two boundary temperatures are only 
indicative when the two temperatures become of order 
T c ] Note also that the region where the AL type con- 
tribution dominates over the BCS bubble is resticted 



to very low temperatures in weak coupling, while it in- 
creases rapidly in the crossover region. In particular, 
it reaches about one-fifth of the critical temperature at 
(kpap)^ 1 = 0, a coupling value for which the chemi- 
cal potential is still positive and remnants of the Fermi 
surface survive (the chemical potential changing sign at 
(kpap)^ 1 = 0.55). The inclusion of the AL type diagram 
turns thus out to be appreciable as soon as one leaves 
the extreme weak-coupling region. In the strong-coupling 
side of the crossover region (i.e., for (kpap)^ 1 > 0.3), the 
AL type contribution definitely dominates over the BCS 
contribution for all relevant temperatures. 




0.1 0.2 0.3 0.4 0.5 



T/T c 

FIG. 4. Temperature dependence of the AL type 
(dashed-dotted line) and BCS (dotted line) contributions to 
the superfluid density, and sum of the two contributions (full 
line), vs temperature when (fc^cr) -1 = 0.0, 0.2, 0.4, from top 
to bottom. [The temperature is in units of T c as obtained 
from Fig. 3.] 

The temperature dependence of the AL type and BCS 
contributions to the superfluid density is shown in Fig. 4 
for three characteristic couplings in the crossover region, 
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namely, (Af^f) 1 = 0.0,0.2,0.4. We have extended this 
calculation up to halfway the critical temperature, where 
the values of A and fi obtained as described above are 
still reliable. Note that for (kpap)- 1 = 0.0 the BCS 
contribution dominates over the AL type contribution, 
while for (Afof) -1 = 0.4 the opposite occurs. This im- 
plies that, even for the supcrfluid density, the crossover 
from weak- to strong-coupling behavior occurs in a rather 
narrow region of the coupling parameter (^pop) . The 
mechanism for this rapid crossover is that the BCS con- 
tribution becomes quickly negligible when (fc^ap) -1 ap- 
proaches zero from negative values ((fcpai?) -1 = cor- 
responds to the threshold for the occurrence of a bound 
state in the associated two-body problem); the AL type 
contribution, on the other hand, becomes quickly negligi- 
ble when (fej?aj?) _1 approaches zero from positive values. 
It turns thus out that there exists a finite interval of the 
parameter (/cfof) -1 where both contributions are small 
with respect to the full density n (at fixed temperature). 
From a physical point of view, one intepretes this result 
by noting that there are two possible mechanisms to de- 
stroy the superfluid density, namely, the pair-breaking 
effect on the weak-coupling side and the excitation of 
collective modes on the strong-coupling side. As both 
mechanisms tend to be suppressed in the intermediate- 
coupling region, the value of the superfluid density is cor- 
respondingly enhanced and superfluidity gets more ro- 
bust against thermal fluctuations. 

To establish comparison with the experimental data on 
cuprate superconductors, the calculation of the tempera- 
ture dependence of the superfluid density should realisti- 
cally include also lattice effects and the d-wave character 
of the gap function. Nevertheless, we expect the above 
physical argument (about the increasing robustness of 
the superfluid density against thermal fluctuations when 
approaching the intermediate-coupling region from both 
sides) to remain valid even in the presence of a d-wave gap 
function. In fact, in the presence of a ci-wave gap func- 
tion the BCS contribution gets suppressed even faster 
for increasing coupling as soon as the chemical potential 
crosses the bottom of the single-particle energy band, 
passing from a linear-T behavior 32 ' 33 to an (s-wave like) 
exponential behavior; the AL type contribution in the 
presence of a d-wave gap, on the other hand, should not 
be appreciably modified with respect to the results ob- 
tained with an s-wave gap, at least when the chemical 
potential remains below the bottom of the single-particle 
band. 

An example of a superconductor with an s-wave gap 
function, for which the effective coupling may be suf- 
ficiently strong that the effect of collective modes on 
the superfluid density could be revealed, is the recently 
discovered MgB^. Experimental data for the temper- 
ature dependence of the penetration depth (and thus of 
the superfluid density) show a quadratic temperature de- 
pendence over most of the temperature range up to T c , 
but for noticeable deviations in the low-temperature re- 
gion which appear instead to follow a quartic behavior 39 . 



While the quadratic temperature dependence away from 
the low-temperature region may simply result from an or- 
dinary (s-wave) BCS contribution (which can indeed be 
rather well fitted to a quadratic behavior upon decreas- 
ing the temperature from T c ), the quartic behavior at low 
temperature may instead be due to the presence of col- 
lective modes. [The occurrence of the above quadratic 
temperature dependence at high enough temperatures 
has alternatively be considered as an indication of the 
presence of nodes in the superconducting energy gap of 
AlgB 2 . 39 } In Fig. 5 we report the superfluid density ver- 
sus T 2 , as obtained by our calculations with both BCS 
and AL type contributions included, for a specially cho- 
sen value of the coupling ((fefaf) -1 = —0.4) such that 
both the T 4 behavior (due to the AL type contribution 
at low temperatures) and the T 2 behavior (due to the 
BCS contribution at higher temperatures) are evident. 
Note that, to evidence the T 2 behavior at high enough 
temperatures, the numerical results reported in Fig. 5 
have been extended beyond their strict range of validity, 
where the values of A and fi should rather be obtained 
by an improved calculation as mentioned above. For this 
reason, the temperature scale in Fig. 5 has been normal- 
ized with respect to the mean-field critical temperature 
Tmf at which A vanishes in the present calculation (and 
not with respect to T c reported in Fig. 3). Note also that, 
while the experimental curve reported in Ref. 39 shows a 
change of convexity when passing from the T 2 to the T 4 
behavior, the theoretical curve shown in Fig. 5 mantains 
the same convexity over the whole temperature range. 




0.1 0.2 0.3 0.4 0.5 

(T/T MF ) 2 

FIG. 5. Superfluid density vs T 2 when (fc f a f ) _1 = -0.4, 
including both BCS and AL type contributions. [The tem- 
perature is in units of the mean-field critical temperature.] 

Care should definitely be exerted when trying to com- 
pare the temperature dependence of the superfluid den- 
sity, obtained also with the inclusion of collective modes, 
with the experimental data for superconductors where 
density fluctuations produce charge unbalance. In this 
case, the Bogoliubov mode with a linear dispersion at 



9 



small wave vectors should turn into a gapped plasmon 
mode owing to the long-range Coulomb interaction. Dis- 
sipation effects have also been recently shown to play an 
important role in restoring a power-law behavior for the 
supcrfluid density at low temperature, albeit with a dif- 
ferent exponent from that obtained by considering only 
neutral collective modes 40 . Renormalization of the chem- 
ical potential (which drives the system from the BCS to 
the BE limit) has, however, not been considered in this 
context. A more realistic approach to the calculation of 
the low-temperature behavior of the superfluid density 
would then include both Coulomb and damping effects 
on top of the treatment presented in this paper, which 
is appropriate to follow the crossover from the BCS to 
the BE limit. Assessment of the relative importance of 
quasi-particle excitations and collective modes to the su- 
pcrfluid density, in fact, can be realistically done only in 
the context of this crossover when the chemical poten- 
tial is allowed to be properly renormalized. This more 
complete approach is deferred to future work. 

Notwithstanding these limitations, it is interesting to 
exploit further our calculation for the current correlation 
at zero frequency, to obtain the form of the Pippard-likc 
kernel from weak to strong coupling at zero temperature. 



B. Pippard-like kernel 

As discussed in the Introduction, two different length 
scales £ p hasc (describing correlation among different 
Cooper pairs or composite bosons) and £ pa ir (describing 
correlation between fermions in a pair) enter the BCS- 
BE crossover problem. In particular, while £ p hasc ~ £ P air 
in the weak-coupling limit, £ p hasc 3> £ P air in the strong- 
coupling limit. 6 In addition, in the weak-coupling limit 
Cpair coincides with the Pippard coherence length £ char- 
acterizing the non-local relation between the supercur- 
rent j and the vector potential A. 26 It is then natural to 
ask which one of the two lengths (£ pa ir or £ p hase) enters 
this relation between the supercurrent and the vector po- 
tential also in the intermediate- and strong-coupling re- 
gions. In this subsection, we consider the contribution 
to the j vs A relation from the standard BCS bubble 
and the AL type diagram introduced in Section III, and 
show that £ pa ; r and £ p h ase are associated with the two 
contributions, in the order. 

Quite generally, to extract the Pippard-like kernel 
(whose precise definition will be given below) from the 
expression of the current correlation function in a form 
appropriate to numerical calculations, it is convenient to 
identify from Eq. (25) the function G AL (q, p), such that 



x ^(Q,O, = 0) = 



(2m) 2 
x (2<2y + Qy) 



(27r)3 (^ + Q,) 

G AL (q,q+Q) 
2|q||q+Q| 



(30) 



for the transverse components (we shall restrict, in par- 
ticular, to the zero-temperature limit). With this defi- 
nition, the function G AL (q, p) reduces to the expression 
(C7) of Appendix C in the strong-coupling limit. 

For all couplings, it turns out that G AL (q, p) = 
G AL {\q\, |p|) and that G AL (q,q) = 0, as in the Bogoli- 
ubov case (cf. Appendix C). [We have verified these prop- 
erties numerically for the AL case within a 5 % accuracy] 
In this way, Eqs. (C9)-(C12) of Appendix C remain for- 
mally valid for all couplings (with the label AL attached 
to the quantities G(p,q) and g(R) therein). 

Strictly speaking, however, it is no longer possible to 
extract the delta-function from the AL contribution at 
arbitrary coupling as in Eq. (C13), since for composite 
bosons \im p ^oo pG AL (q,p) = (instead of the behavior 
(C15) for point-like bosons). For composite bosons the 
coefficient X of Eq. (C13) vanishes and only the nonlo- 
cal contribution g AL (R) survives. Nonetheless, upon ap- 
proaching the extreme strong-coupling limit (where the 
composite nature of the bosons no longer matters) the 
small- R behavior of g AL (R) can be effectively assimilated 
to a delta- function (since its range extends over a length 
scale much smaller than the scale £ p hasc associated with 
the remaining part of g AL (R)), such that the expected 
behavior for point-like bosons is properly recovered in 
this limit, as also shown below by the numerical results. 

Similarly, the BCS bubble contribution can be calcu- 
lated at arbitrary couplings according to the procedure 
of Appendix D, instead of relying on the standard BCS 
analytic calculation for weak-coupling 2 . In this case, one 
can identify the function <?„ s (i?) for R ^ and its in- 
tegral X BCS = -87rW (cf. Eq. (Dll)). What happens 
to the BCS bubble contribution is that, in the weak- 
coupling limit, the function g® cs (R) contains a short- 
range part extending over the scale fc^ 1 <C £ pa ir, whose 
integral (from R = to R = kp 1 ) together with the delta- 
function contribution proportional to X BCS cancels the 
diamagnctic contribution to the j vs A relation, making 
only the nonlocal part of g^ cs (R) effectively surviving. 

Putting together the results for the AL and BCS con- 
tributions, the desired relation between the induced cur- 
rent j and the vector potential A acquires the form: 



n 



+ 



3mc 
1 



R^EyRy'Ay^) 



(27r) 4 m 2 c J R 4 
x(g* cs (R)+g AL (R)) 



(31) 



where R — r r' and c is the light velocity. This re- 
lation defines the Pippard-like kernel at arbitrary cou- 
plings, which is the sum of a local (London-like) and of a 
nonlocal component. Note that the local part of Eq. (31) 
coincides with the local part of Eq. (D12) since the AL 
contribution lacks the delta function. 

The results of the numerical calculation for the two 
(BCS and AL) contributions to Eq. (31) show that, in 
the weak-coupling limit, the BCS contribution dominates 
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over the AL counterpart even though both contributions 
extend over the same spatial range since £ pa i r ~ Cphasc 
in this limit. In the strong-coupling limit, on the other 
hand, the BCS bubble no longer contributes significantly 
to the response function since the nonlocal part g^ cs (R) 
cancels the delta-function contribution over a scale £ pa ; r 
much smaller than the scale £ p hase of the AL contribution. 
At arbitrary coupling, both BCS and AL contributions 
to Eq. (31) need be explicitly considered. 
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FIG. 6. Function g£ L (R) (in units of m B /(4£* haso )) vs R 
(in units of %/2£ p hase) f° r four different values of the coupling 
parameter (&fOf) _1 : —1.0 (dotted line), 0.0 (long-dashed 
line), +1.0 (dashed-dottcd line), +2.0 (full line). The purely 
bosonic g% (R) (short-dashed line) is also shown for compari- 
son. 

In Fig. 6 the functions g AL (R) is reported vs R = 
i?/(\/2Cphaso) for different values of the coupling pa- 
rameter (fop op) -1 . The relevant range of g AL (R) ex- 
tracted from this plot is determined by the coher- 
ence length £ p hasc introduced in Ref. 6 to account for 
the spatial correlations of the amplitude of the or- 
der parameter (direct calculation 37 of fc,F£ P hasc yields 
the values (1.7,0.7,0.8,1.1) corresponding to the val- 
ues (kpap)- 1 = (-1.0, 0.0, +1.0, +2.0) here considered). 
[We have also verified that /cpCphaso extracted from the 
behavior of g AL (R) attains its minimum at (&fof) _1 = 
+0.4, in agreement with the results from Ref. 6.] Note 
that the delta-function contribution to g(R), which is 
present in the bosonic expression of Appendix C, is re- 
covered in Fig. 6 by the progressive shrinking of the short 
range part of g£ L (R) as the coupling is increased; at the 
same time, the long-range part of g AL (R) in Fig. 6 ap- 
proaches the asymptotic behavior of the bosonic g^(R). 
We should warn, however, that extracting the short- 
range part of g^ L (R) from numerical calculation of the 
AL type diagram is affected by sizable numerical errors, 
since it requires an accurate treatment of the large wave- 
vector behavior of G AL (q,p) (this procedure, in turn, re- 
quires subtraction of large terms which should add up to 
zero). 
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FIG. 7. Function g% cs (R) (in units of mk%) vs R (in 
units of £ P air) for four different values of the coupling pa- 
rameter (fefap)" : —2.0 (dotted line), —1.0 (dashed line), 
0.0 (dashed-dotted line), +1.0 (full line). 

In Fig. 7 the function g^ cs (R) is reported vs R = 
-R/£pair for different values of the coupling parameter 
(kpap)^ 1 . The relevant range of g^ cs (R) is deter- 
mined by the coherence length £ pa ir for two-fermion 
correlation 5,6 , which is monotonically decreasing for in- 
creasing coupling. In particular, direct calculation 37 
of fcp^ pa i r yields the values (15.2,3.4.1.1,0.6) for 
(kpapy 1 = (-2.0, -1.0, 0.0, +1.0), respectively, in 
agreement with the values extracted from Fig. 7. 

We recall at this point that (apart from a numerical 
factor of order unity) £ p hase and £ pa i r coincide with each 
other in the intermediate-to-weak-coupling side of the 
crossover (which is mostly relevant to the BCS contri- 
bution). We thus conclude from the present analysis of 
the current correlation function that it is the length £ p hase 
to determine the range of the Pippard-like kernel over the 
whole crossover range. In addition, the present discus- 
sion shows the importance of including both BCS and 
AL type contributions in the intermediate-coupling re- 
gion, for a correct description of the current correlation 
function. 



V. CONCLUDING REMARKS 

In this paper, we have calculated the current corre- 
lation function within the BCS-BE crossover. We have 
complemented the standard analysis of the BCS diagram- 
matic contribution to the current correlation function 
in the broken-symmetry phase, by considering an addi- 
tional (Aslamazov-Larkin type) diagrammatic contribu- 
tion which includes fluctuation effects in weak coupling 
and recovers the Bogoliubov form of the current correla- 
tion function in strong coupling. 

In view of application to high-temperature cuprate su- 
perconductors, the analysis presented in this paper is 



11 



necessarily a preliminary one, since it relies on a contin- 
uum model in three dimensions and does not take into 
account additional effects (such as the Coulomb repul- 
sion between charged composite bosons), which modify 
the Bogoliubov picture for bosons. In addition, lattice ef- 
fects and the <i-wave character of the gap function should 
also be included in a realistic calculation of the superfluid 
density. 

The analysis presented in this paper is, nonetheless, a 
required preliminary step toward a more realistic the- 
ory of the current correlation function in the super- 
conducting phase, by connecting results valid in the 
weak- and strong-coupling limits in the spirit of the 
BCS-BE crossover. In this context, the identification 
of the Aslamazov-Larkin type diagram(s) in the broken- 
symmetry phase as the dominant one(s) in the strong- 
coupling limit of the theory represents one of the main 
contributions of the present paper. In addition, con- 
sideration of the evolution of the current correlation 
function in the broken-symmetry phase from weak to 
strong coupling has required us to deal with several rel- 
evant theoretical issues, such as the mapping between 
the fermionic and bosonic diagrammatic structures in the 
broken-symmctry phase, the connection between the Bo- 
goliubov propagators for composite bosons and the di- 
agrammatic structure for the constituent fermions, and 
the identification of the local and nonlocal parts of the j 
vs A relation for a generic diagrammatic contribution to 
the current correlation function. The temperature depen- 
dence of the superfluid density that we have obtained (by 
including the AL type diagram on top of the conventional 
BCS "bubble" diagram with an s-wave gap function) has 
served to identify the temperature regions where one of 
the two contributions dominates over the other one in dif- 
ferent coupling regimes. In this context, we have shown 
that, in the intermediate-coupling region, superfluidity 
is most robust against suppressing mechanisms, such as 
pair breaking (coming from the weak-coupling side) and 
collective excitations (coming from the strong-coupling 
side), since the effects of these mechanisms is consid- 
erably reduced in the intermediate-coupling region. It 
would accordingly be tempting to relate this effect to the 
occurrence of a maximum for the superconducting criti- 
cal temperature in the phase diagram of cuprate super- 
conductors at intermediate doping. We have also been 
able to extract from our (zero-temperature) calculation 
of the spatial dependence of the current response kernel 
a length scale previously introduced in the theory of the 
BCS-BE crossover, which specifies the correlation among 
different fcrmion pairs. 

Although detailed comparison with experimental data 
on cuprate superconductors has to await more realistic 
calculations that include lattice effects and the ci-wave 
character of the gap function (for instance, via an ex- 
tended negative-C/ Hubbard model), the trend of the co- 
herence length in the Pippard-like kernel to decrease ini- 
tially with increasing coupling from weak to intermediate 
should be experimentally accessible for cuprate supercon- 



ductors upon decreasing the doping (from ovcrdoping to 
underdoping) , in analogy with the original experimental 
results by Pippard 26 . 

Previous calculations of the temperature dependence 
of the superfluid density have included the ci-wave charac- 
ter of the gap function within the standard BCS "bubble" 
contribution 33 , or have considered incoherent pair exci- 
tations of finite momentum which assist the destruction 
of superconductivity (still taking into account the ci-wave 
gap) in the context of the BCS-BE crossover 41 . As both 
works do not address the role of collective modes (over 
and above the quasi-particle contribution) in the context 
of this crossover, direct comparison with our work cannot 
be made. Additional inclusion of long-range Coulomb in- 
teraction addressed in Ref. 40 is more directly relevant to 
the present approach and is presently being considered. 

ACKNOWLEDGMENTS 

The authors are indebted to C. Castellani for discus- 
sions. Financial support from the Italian MIUR under 
contract COFIN 2001, Prot. 2001023848 is gratefully ac- 
knowledged. 



APPENDIX A: PARTICLE-PARTICLE LADDER 
PROPAGATOR IN THE BROKEN-SYMMETRY 
PHASE AND THE BOGOLIUBOV 
PROPAGATORS 

In this Appendix, we study in detail the particle- 
particle ladder entering the fermionic diagrammatic 
structure for the broken-symmetry phase, and specifi- 
cally the AL diagram(s) considered in Section III. To 
this end, we shall take the fermionic attractive interac- 
tion of the form of a contact-potential with a suitable 
rcgularization, as it was done in Ref. 9 for the normal 
phase. In addition, in the strong-coupling limit we es- 
tablish a connection between the particle-particle ladder 
propagator and the Bogoliubov propagators for bosons, 
correcting at the same time some erroneous conclusions 
reached in Ref. 4 in this context. 

Knowledge of the detailed form of the attractive in- 
teraction V(r) of Eq. (12) is not required when studying 
the BCS-BE crossover. One may accordingly consider 
the simple form of a "contact" potential V(r) = voS(r), 
where vo is a negative constant. This choice of the po- 
tential entails a suitable regularization, e.g., in terms of a 
cutoff fc in wave-vector space. In three dimensions, this 
can be achieved via the scattering length aj? of the as- 
sociated two-body (fermionic) problem, by choosing the 
constant vq as follows 9 : 



mk mapkQ 
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In this way, the classification of the (fermionic) many- 
body diagrams is considerably simplified even in the 
broken-symmetry phase, since only specific diagram- 
matic substructures survive when the limit k — > oo is 
eventually taken. 



and 



p+q 



P +q p+q p'+q p+q p'+q 



r(q) 



1 


2 

2 2 1 






r 22 







■• (a) 



(b) 



(c) 



FIG. 8. (a) Particle-particle ladder in the 
broken-symmetry phase (conventions for four-momenta and 
Nambu indices are specified); (b) and (c) show typical inser- 
tions of the rungs II12 and II21, which are compatible with 
the regularization (Al) of the fermionic interaction potential. 
[Dots delimiting the potential (broken line) here represent tz 
Pauli matrices.] 

With these premises, we consider the particle-particle 
ladder depicted in Fig. 8(a). 30 One can convince oneself 
that, in these diagrams, the elementary rung 



Ki L e'£ R e'(q) 



■ q)Ge> R e> L (-p) 



(A2) 



with the single-particle Green's functions of the BCS 
form (13), survives the regularization introduced by the 
potential (Al) in the limit k n — > 00 (and thus v — ► 0) 
only if£ R = £ L =£' R =£' L = 1 or £ R = £ L = £' R = £' L = 2. 
[We use the coventions of Fig. 1 for the Nambu indices.] 
We then isolate the ultraviolet divergence of (A2) as 
k — > 00, by writing [cf. Eq. (17)] 



Hii (q) 
as well as 



# 11,11 (9) = Ru(q) + 



mho 
2^" 



1122(g) = K 2 i,2i{q) = nu(-g) . 



(A3) 



(A4) 



where 



Rn{q) 



m 

n(-P)-0 



is a well-behaved quantity in the limit fco — > 00. 
way, we are led to introduce the partial ladders 

f n(g) = -v + ^o n u(g) - «o n ii(g) 2 + • • 
1 



4ttcif 



+ Ru{q) 



k 2 

(A5) 
In this 



(A6) 



f 22(g) = -v + vln 22 (q) - -yon 22 (g) 2 + • 
= f n (-g). 



(A7) 
'a) there 



£' L different 



This does not imply, however, that in Fig 
cannot occur structures with the pair £l = 
from the pair £r — £' R . The rungs ni 2 (g) = ^11, 22(g) or 
n 2 i(g) = ^22, 11(g) [cf. Eq. (18)] can, in fact, be inserted 
in between the above partial resummations r(g)n and 
T(g)22, as indicated in Figs. 8(b) and 8(c). It is possible 
to take into account these structures in a systematic way, 
by introducing the matrix equation 



Tn(g) r i2 (g) 
T 2 i(g) r 22 (g) 



-w 



nn(g) 
n 2 i(g) 



n i2 (g) 
n 22 (g) 



1 
1 

rn(g) 
r 2 i(g) 



(A8) 



Fi 2 (g) 
r 22 (g) 



whose solution yields 



Tn(g) T 12 {q) 

r 2 i(g) r 22 (g) 



A(q)A(-q) - B{qY 

M-q) B(q) \ 

B{q) A(q) ) 



(A9) 



with the short-hand notation 
-A(q) = - + n n (g) 

B(q)=U 12 (q). 
Note that in Eq. (A10) the gap equation 

dp 1 



v J (2ir 



(2tt)3 2E(p) 



2E(p) 
(A10) 



(All) 



has been taken into account. In this way, we can repro- 
duce all diagrams in Fig. 8 that survive the regularization 
of the potential. 

The expressions (A10) for A(q) and B(q) simplify con- 
siderably in the strong- coupling (/3/i — > —00) limit we are 
interested in. One obtains 37 : 



A{q) a 



Ait 



apkp 
3tt 



and 



B{q) 



m 



L F F 



q 2 - 
4 8 



12^ 2 



(A12) 



(A13) 



to the lowest significant order (where the Fermi wave 
vector fop is related to the density n via n = k F /(3n 2 )). 
This gives for the matrix elements (A9) the expressions: 
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rn(g) = r M (- g) * + ;"; 2 + ? V( l m) (aw) 



m 2 as -Es(q) 2 - (iw„) 2 



and 



ri 2 (g) = r 2 i( g )~ 



8tt 



Mb 



where 



#u(q) = 




Ms -Ms 



(A15) 



(A16) 



nig = 2m being the bosonic mass and Ms = 
2k F a F I (3nm) the bosonic chemical potential. Note that 
Ms can be cast in the Bogoliubov form 



Ms = n B v 2 (0) = 



n 4ira F 



m 



(A17) 



in terms of the bosonic density = n/2 and of the 
residual bosonic interaction 112(C)) = Ana F /m. ' 9 The ex- 
pression (A16) has thus the form of the Bogoliubov dis- 
persion relation 2 . 

The results (A14) and (A15) bear strong resemblance 
with the normal (£?') and anomalous (G'21) nonconden- 
sate bosonic Green's functions, respectively, within the 
Bogoliubov approximation, which read 2 : 



<?'(<?) 



Ms + iu v + q 2 /(2ms) 



(iu v ) 2 - E B (q) 2 



Ms 



(iw„) 2 -£fl(<0 



2 • 



(Ais) 

(A19) 



Comparison of Eqs.(Al4)-(A15) with Eqs. (A18)-(A19) 
thus lead us, in fact, to identify: 



r n (g) = -(8n/m 2 a F )G'(q) 
Ti 2 {q) = (8n /m 2 a F )g' 21 (q) 



(A20) 



where the sign difference between the two expressions 
must be remarked. 

Haussmann 4 had already noted this sign difference and 
attempted to get rid of it by fixing the phase of the 
BCS gap function at the particular value tt/2. This pre- 
scription appears rather arbitrary, since the theory can- 
not depend on the value of the phase of the gap func- 
tion. This sign difference is instead genuine, reflecting 
the fact that the particle-particle ladder propagator does 
not fully reduce to the Bogoliubov propagators in the 
strong-coupling limit. 

To demonstrate this point explicitly, we resort directly 
to the definition of the Bogoliubov propagators 2 

G'(x,x') = -(T T [* B (x)* B (x')])+a(r)a*(r') (A21) 

and 

g' 21 (x,x') = -{T r [*Ux)*Ux')}) +a*(r)a*(r% (A22) 



where x = (r,r) and a(r) = (^s(r)) is the condensate 
amplitude defined in terms of the bosonic field operator 
^s(r). The desired connection between fermionic and 
bosonic quantities then results by relating the bosonic 
operator \1/ B (r) to its fermionic counterpart Vv (r) (where 
°~ = (Tj I) i s the spin projection), by setting 6 : 

tf B (r) = J dpcb(p)^ L (r - p/2)^ T (r + p/2) (A23) 

where 4>(p) is a suitable function. On physical grounds, 
we choose <j>(p) to be the (normalized) bound-state wave 
function of the associated (fermionic) two-body problem: 



-p/a F 



^2ita F p 
with Fourier transform 

<Hp) = 



a F p 2 + a F 2 



(A24) 



(A25) 



Entering the definition (A23) into Eqs. (A21) and (A22), 
and recalling the definition (5) of the two-particle corre- 
lation function, we rewrite the Bogoliubov propagators 
in the compact form: 

g'(rr,r'r') = - Jdp j dp' <j>(p)<}>* (p')L(l, 2; l', 2') (A26) 

with the identification 1 = (r + p/2,r,£ = 1), 2 = (r' — 
p'/2, t', I = 2), 1' = (r - p/2, r+, I = 2), and 2' = (r' + 
P '/2,t'+,£= 1); as well as 

g' 21 (rr, r'r') - - J "dp J dp' ^ (p)<f>* (p')L(l, 2; l', 2') 

(A27) 

with the identification 1 = (r — p/2,r,£ — 2), 2 = (r' — 
p'/2, t', I = 2), 1' = (r + p/2, t+,£=1), and 2' = (r' + 
p'/2,T ,+ ,£ =1). In Fourier space, we write: 

x <t>(p' + q/2)L^(pw„,pW; q) (A28) 

and 

x 0(p' + d/2)L 21 2 (pu> pW;g) (A29) 

where the conventions for the four-momenta and Nambu 
indices are as in Fig. 1. 

We can use at this point the expression (9) for L in 
terms of the many-particle T-matrix, and limit ourselves 
to consider the BCS approximation (12) for the kernel S, 
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as obtained from the off-diagonal terms of the BCS self- 
energy. As already stated in Section II, with this restric- 
tion only four elements of the many-particle T-matrix 
with £l 7^ £'l and £r ^ £' R survive. In addition, for a 
contact potential the T-matrix in Fourier space satisfies 
Eq. (16), whose solution yields 



Tu(q) T 12 (q) 
T 21 (q) T 22 {q) 



A{q)A{-q) - B(q) 2 

-A(-q) B{q) 
B{q) -A(q) 



(A30) 



with the same definitions of Eq. (A10). 

We are specifically interested in the strong-coupling 
limit of the above expressions, whereby the off-diagonal 
single-particle Green's functions can be neglected in com- 
parison to the diagonal ones since A <C \p\ in this limit. 
In this way, the expressions (A28) and (A29) reduce to: 



g'(q)^-T 2 {q)+T 1 (q) 2 T 11 {q) 



(A31) 



and 



021 (9)^1 (-9)T2i(g) (A32) 

where we have introduced the notation (j = 1,2) 

^ {q) = j ^<Kp + q/2)^E^ + ^n(-p)- 

(A33) 

With the approximate result 
jj$^0ii(p + «)0ii(-p) 



£(k) + E(k + q) - %u v 
1 



2E(k) 



(A34) 



that holds specifically in the strong coupling- limit, and 
using the expression (A25) for the two-body wave func- 
tion, Eq. (A33) can be readily evaluated to yield: 



Ml) 



m 2 ap ^_ / % ma 



(A35) 



Neglecting the first term on the right-hand side of 
Eq. (A31) (which does not contain the polar structure 
of the T-matrix), we obtain eventually: 



, m 2 a F 
5' 21 (q)- -^TM- 



(A36) 



This result show that (apart from a constant prefactor) 
it is the T-matrix and not the particle-particle ladder 
to reduce to the Bogoliubov propagators in the strong- 
coupling limit of the fcrmionic theory. In particular, com- 
parison of Eq. (A9) with Eq. (A30) yields Eqs. (19) of 



the text, thus accounting for the sign difference noted in 
Eqs. (A20). 

Finally, it is interesting to evaluate in the strong- 
coupling limit also the condensate amplitude entering 
Eqs. (A21) and (A22), with the same approximations 
used above. From its definition in terms of the bosonic 
operator (A23), we obtain (for a homogeneous system): 



a = jdp^p)^^- p/2)^{v + p/2)) 

A r dp 0( P ) 

2 J (2tt) 



r)3 E(p) 



m 2 a,F 
8tt 



(A37) 



Note that a similar relation holds (with a in Eq. (A37) re- 
placed by y/n/2) when the diagonal BCS single-particle 
Green's function Gn of Eq. (13) is used to evaluate the 
particle density n in strong coupling and at low temper- 
ature 37 . 



APPENDIX B: MAPPING OF THE FERMIONIC 
ONTO THE BOSONIC DIAGRAMMATIC 
STRUCTURE IN THE BROKEN-SYMMETRY 
PHASE 

In Appendix A we have established a connection 
[cf. Eqs. (A20)] between the Bogoliubov propaga- 
tors (A18)-(A19) and the particle-particle ladder (A14)- 
(A15), in the strong-coupling limit. Since the Bogoli- 
ubov results (A18)-(A19) for true bosons can be obtained 
by considering self-energy corrections to the free-boson 
propagator 35 ' 36 ' 2 , it should also be possible to obtain the 
expressions (A14) and (A15) for the particle-particle lad- 
der in the strong-coupling limit, by considering bosonic 
self-energy corrections (expressed in terms of the effective 
bosonic interaction in the normal phase and of the con- 
densate density) to the particle-particle ladder in the nor- 
mal phase. In this way, a mapping between the fermionic 
and bosonic diagrammatic structures is effectively estab- 
lished in the broken-symmetry phase, in a similar fashion 
to what has already been achieved in the normal phase 9 . 

We begin by recalling that the Bogoliubov propa- 
gators (A18)-(A19) in the broken-symmetry phase can 
be obtained from the free-boson propagator Q^°\q) = 
(iujv — q 2 /(2ms) 2 + ps) 1 , by considering the following 
normal and anomalous bosonic self-energy corrections: 



Sf 1 (g) = Sf 2 (- 
Ef 2 (g) = S 2 B 1 (- 



q)=n o (V B (0) + V B (q}) 
q) = n V B (q) 



(Bl) 



where Vs(q) is the Fourier transform of the bosonic in- 
terparticle potential. Within the same approximation, 
the (bosonic) chemical potential is given by 



p B = ££(0) 



Sf 2 (0) 



n o V B (0), 



(B2) 
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in agreement with the Hugcnholtz-Pines theorem 2 . 

To obtain the corresponding self-energy corrections for 
composite bosons, it is sufficient to consider the rungs 
Ilii(g) [cf. Eq. (17)] and U 12 (q) [cf. Eq. (18)] in the 
strong-coupling limit. At the lowest significant order, we 
obtain for the contributions proportional to A 2 : 



n ii(<z) — > -n ma 2 F 

ni 2 (<?) ^ Y maF 



(B3) 



where n = \a\ 2 is the condensate density. Recalling that 
in the normal phase the particle-particle ladder reduces 
to 4 ' 9 



r°(<f) - - 



i 



8tt 



m 2 aF ) iw v — q 2 /(4m) 2 + /is 



(B4) 



in the strong-coupling limit, the expression 
taining a single insertion of IIn reads: 



for Tn con- 



r°(g)n u ( g )r ( ? ) 

8tt 



m 2 ciF J iu v — q 2 /(4m) 2 + [Ib 
8tt \ 1 



(-n ma 2 F ) 



"m?ap ) iuj v — q 2 /(4m) 2 + hb 
8tt \ 1 



m 2 ciF J iu v — q 2 /(4m) 2 + hb 
1 



(2n v 2 (0)) 



iuj v - q 2 /(4m) 2 + [i B 



(B5) 



where v 2 (0) — AimF/rn is the residual bosonic 
interaction 4 ' 9 introduced in Eq. (A17). Apart from the 
overall factor — Sir / (rn^ap), the expression (B5) for com- 
posite bosons coincides with the lowest-order diagram 
obtained within the Bogoliubov approximation for true 
bosons, with the self-energy insertion £n(g) = 2n v 2 (Q) 
(which coincides with the first of Eqs. (Bl), provided one 
approximates Vs(q) ~ Vb(0) — > v 2 (Q) for the relevant 
range of wave vectors). 

By a similar token, consideration of the structure of 
r (g)rii2(<7)r (-g) leads us to identify Y, 12 {q) = n v 2 (0) 
as the off-diagonal self-energy insertion for composite 
bosons to lowest order in A 2 . [Note that it is just the 
sign difference between Ilnfg) and Ui 2 (q) in Eqs. (B3) 
which eventually results in the sign difference between 
Tn and T\ 2 with respect the corresponding Bogoliubov 
propagators, as evidenced in Eqs. (A20).] 
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<=> X/ -k + q r q 4 



(a) 
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FIG. 9. Correspondence between the: (a) Effective 
two-boson interaction and the four-point fermionic interac- 
tion (with qi +q2 = q3 + q<i); (b) Anomalous and (c) normal 
self-energies within the Bogoliubov approximation and the 
contractions of the four-point fermionic interaction of (a). 

The above results can be interpreted euristically in 
terms of an effective interaction between composite 
bosons in the broken-symmetry phase, as follows. Let's 
consider the four-point vertex for composite bosons in the 
normal phase depicted in Fig. 9(a), which plays the role 
of the effective boson-boson interaction 4 ' 9 in the strong- 
coupling limit (we shall consider a definite choice of the 
spin indices). 

To obtain the bosonic self-energy Sf 2 of Fig. 9(b) one 
should: (i) contract the two pairs of "outgoing" connec- 
tions in Fig. 9(a), by joining them through multiplication 
with the BCS value of —A; (ii) multiply each fermionic 
single-particle line proceeding in the opposite direction 
(with respect to one pair of "ingoing" connections) by 
— 1, thus matching the standard convention [cf. the first 
of Eqs. (13)]. The final result corresponds to ITi2(^) of 
Eq. (B3), since (to the lowest order in A/|/x|) Q\ 2 can be 
written in the form: 



Qii{q) ^ 



(-A)- 



(-1) 



■Ju - £(k) 



(B6) 



Note that in the above equation A plays the role of the 
square root of the condensate density (apart from a nu- 
merical coefficient), in agreement with Eq. (A37). 

By a similar token, to obtain the bosonic self-energy 
Sf x of Fig. 9(c) one should: (i) contract one pair of "out- 
going" and one pair of "ingoing" connections in Fig. 9(a), 
by joining them through multiplication with the BCS 
value of —A in two possible ways; (ii) multiply again 
each fermionic single-particle line proceeding in the op- 
posite direction (with respect to one pair of "ingoing" 
connections) by —1. In the strong-coupling limit, the 
final result corresponds to — 2IIii(g = 0) of Eq. (B3). 
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APPENDIX C: EXTRACTING THE PIPPARD 
KERNEL FOR BOSONS WITHIN THE 
BOGOLIUBOV APPROXIMATION 

In this Appendix, we obtain the (wave-vector and 
frequency dependent) bosonic current correlation func- 
tion within the Bogoliubov approximation, to which 
the asymptotic form of the fcrmionic current correlation 
function in the strong-coupling limit needs be compared. 
We shall obtain that by relying on a diagrammatic de- 
coupling method in the broken-symmetry phase 42 . The 
static expression of this correlation function is then an- 
alyzed to extract the nonlocal part of the j vs A rela- 
tion (Pippard kernel) . Finally, the penetration depth ob- 
tained from Bogoliubov theory is briefly discussed. 

In analogy to Eq. (10) for fermions, the following ex- 
pression holds for the bosonic current correlation func- 
tion: 



X* y (rr,rV) 



1 



{2m B } 



f (V 7 -V 7 )(V 7 ,-V 7 ,) (CI) 



x (T T [*| J (rT+)* i j(r"T) 

X *^rV' + )* S (r'"T')])|r"=r,r-=r' 

in terms of the bosonic field operator ^ B (r) introduced 
in Appendix A. Setting further 



* B (r) = a(r) + ip{r) 



(C2) 



where a(r) = (^ , s(r)) is the condensate amplitude and 
ip (r) the deviation operator 2 , the right-hand side of 
Eq. (CI) is equivalent to the sum of 16 time-ordered 
products. With the factorization 

(T r b(lM2Vt ( lV(2')]) 
H^(lM2)])(T T [^(iV(2')]) 
+<T T [ V (l) V t(i')])(T T [ V (2) V t( 2 ')]) 
+<T T [ V (l) V t( 2 ')])(T T [ V (2) V t(i')]) (C3) 

which holds for a system of non-interacting (bosonic) 
quasi-particles (as described by the Bogoliubov approxi- 
mation) , and restricting to the case of a uniform conden- 
sate whereby a(r) = a, the Fourier transform Eq. (CI) 
becomes eventually: 



A 7 , 7 



(Q,fi„) = f d(r-T')e in "< T - T '> f d(r - v') 
xe- iQ ' (r - r,) X 7 S f (rr,rV) 



i 



a 



"' 7 (2m B ) 2 
dq 



[G'(Q) + g'(-Q)-2g' 21 (Q)} 

j(2g 7 + Q 7 )(2 9y +Qy)4X! 



(2m B ) 2 J (2tt)3^ ' '"VP 

x[g'(q)g , (q + Q)-g , 21 (q)g , 21 (q + Q)]. 



(C4) 



Utilizing at this point the form (A18 ) and (A19) of the 
Bogoliubov propagators, the result (C4) is seen to co- 
incide with the one obtained alternatively by exploiting 



the Lehmann representation of the current correlation 
function 42 . Although corrections to the Bogoliubov ap- 
proximation for the current correlation function (as to 
make the longitudinal f-sum rule satisfied while preserv- 
ing the form of the Bogoliubov propagators) have been 
considered in a formal way 43 , in this paper we restrict 
ourselves realistically to the simplest approximation (C4) 
for the bosonic limit, since we are ultimately interested in 
the evolution of the fermionic current correlation func- 
tion from weak to strong coupling. 

The "static" (Q -> 0, Q, v = 0) limit of Eq. (C4) can be 
used to test explicitly the longitudinal f-sum rule (15) and 
to obtain the (temperature dependence of the) superfluid 
density via Eq.(14). One obtains: 



x 7>7 '(Q 



o,n„ = o) 

n 



lim QyQ 



(u B (Q)+v B (Q)f 



+ 



2m|Q^o^ 7 ^ 7 E B (Q) 
dq g 7 q y dn B (E B (q)) 
(2n) 3 m B m B dE B (q) 



(C5) 



where n B (E) = {e@ E — 1) 1 is the Bose-Einstein 
distribution function, 2v B (q) 2 = [(q 2 /(2m B ) + 
n V B (q))/E B (q) - 1] and u B (q) 2 - v B (q) 2 = 1, V B (q) 
here being the same bosonic potential of Appendix B. 
Recalling that v B (q) 2 — ms/(2\q\) and E B (q) ~ s|q| in 
the small- q limit (where s = \/noV B (0)/m B is the sound 
velocity) , the first term on the right-hand side of Eq. (C5) 
reduces to — Q 7 Q 7 'no/raB, implying that the f-sum rule 
(15) is only approximately verified, insofar as no ~ n B 
within the Bogoliubov approximation. The second term 
on the right-hand side of Eq. (C5), on the other hand, 
has the form of {—l/m B times) the expression of the Lan- 
dau's normal- fluid density p n = n — p s (cf., e.g., Ref. 2), 
and is also temperature dependent (albeit through the 
Bogoliubov expression of the quasi-particle excitation en- 
ergy). 

In addition, the Pippard kernel can be obtained from 
Eq. (C4) by letting — but keeping Q finite, in anal- 
ogy to the BCS treatment for fermions 1 ' 2 . In particular, 
for a transverse vector potential in the zero-temperature 
limit one considers the expression: 

J(#H R - Q x 7 VQ^ = o) 



2m% 



j p OO p OO 

——j / dqq / dppG B (q,p) 



r^o dw: joiq]R ' R|) - ?o(p|R ' + R|) (C6) 



7 

with the notation q = \q\,p = |p|, and 

qp 



G B {q lP ) = 



ltB(q)eB(j>) 



E B (q)E B (p)(E B (q) + E B {p)) 1 
+ n V B {q)e B (p) +n V B (p)e B (q) - E B (q)E B (p)} (07) 

where e B (q) = q 2 /(2m B ) and jo(z) — (smz)/z is the 
spherical Bessel function of zeroth order. The function 
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G B (q,p) (plotted in Fig. 10) is symmetric under the in- 
terchange of q and p, vanishes for q — p but is elsewhere 
positive, and its relevant range extends about q ~ Cphasc 
and p • 



€ p iLe> where 



£pha 



1 



(C8) 



is the characteristic length of the Bose gas associated 
with the chemical potential [1b = «oVb(0) [cf. Eq. (A17)] 
which represents the only energy scale in the problem 44 . 




FIG. 10. Two-dimensional plot of the function G B (q,p) 
defined by Eq. (C7), with the wave vectors in units of 

se J 



For the purpose of extracting the leading behavior of 
the function (C6) for large values of R = |R|, the spatial 
derivatives therein can be calculated with the same ap- 
proximations exploited in the analogous calculation for 
the BCS (bubble) contribution to the current correlation 
function 1,2 . One obtains: 



lim 

R'^0 



^^rjo(«|R'-R|)jo(p|R' + R| 
-2^cos(q- P )R. 



(C9) 



Introducing the notation 

pOO pOO 

g B (R)= dqq dppG B (q,p)cos(q-p)R, (CIO) 



Jo Jo 
Eq. (C6) acquires eventually the form: 

/r§ eiR ' Q ^' (Q '^ = 0) 



R>yR. 



(27r) 4 m|j i? 4 



g B {R). 



(Cll) 



It should be noted that 

fOO 



poo poo 

/ dRg B (R) = ir dqq 2 G B {q 1 q) = (C12) 
Jo Jo 

since G B (q, q) = for all values of q. 

In practice, it is convenient to extract at the ouset a 
delta-function contribution from the function g B (R), by 
setting 



where 



g B (R) = g B (R) - 21 B S(R) 



1 B = I dRgZ(R). 
lo 



(C13) 



(C14) 



and evaluate the remainder g B (R) numerically. To this 
end, we exploit the asymptotic form of the function 
G B (q,p) and set 



f B (q) = Jim pG B {q,p) = 2m B 

p—>oo 



(V + 1 - V¥+W) 



7,2 



(C15) 



where q = gV2£ ph asc and p = pV2£ ph a S c- With the fur- 
ther notation 



G B (q,p) = qG B (q,p)p-qf B (q)-pf B (p), (C16) 



z tBf~\ 



we obtain: 



g B n {R) 



l 



4e 4 



dq dpG (q,p) cos(q — p)R 

phase wO J 



2_ r 

RJo 



f-/ dqqf B (q)sm(qR) 



where R — ii/(\/2£phase)> as well as 



irV2 

^•phasc 



pOO 

/ dqqf B (q) 
JO 



nrriB 



se 3 



(C17) 



(C18) 



phase 



The function g B (R) (plotted in Fig. 11) is negative, di- 
verges like tob(7t— l)/(2£ p haso) In .R for R — > 0, is strongly 
localized over the range £ p hase, and shows minor oscilla- 
tions about zero for R > £ p hasc- 

In contrast to an analogous procedure within the 
fermionic BCS approximation of Appendix D, the 
strength X B of the delta function in Eq. (C13) depends 
on the condensate density via £ p hase (in particular, the 
strength of the delta function and consequently the area 
enclosed by g B vanish in the limit £ p hasc — » oo, cor- 
responding to noninteracting (condensed) bosons). Ac- 
cordingly, the delta-function contribution does not cancel 
the diamagnetic term in the relation between the induced 
current j and the vector potential A, as it happens in- 
stead for the BCS case 2 (cf. also Appendix D). In the 
present case, in fact, one obtains: 
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i 7 ( r ) 



1 



m B c 
1 



n B 



3fi7T 2 £ 3 
OU7r ?phasc 



A 1 (r) (C19) 



(27r) 4 m|c 



dr 



i? 4 



where R = r — r' and c is the light velocity. Note how the 
form (C19) for the current response kernel is intermediate 
between the (local) London's form and the (nonlocal) 
Pippard's form 2 . Note also that, in the noninteracting- 
boson limit (whereby £ p hase — » oo), Eq.(C19) correctly 
reduces to the local London's form 



(C20) 



in agreement with the result obtained in Ref. 45 for a 
condensed Bose gas at zero temperature (for fermions 
within the BCS approximation, on the other hand, only 
the nonlocal Pippard's term survives 2 ). 
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FIG. 11. Plot of the function gS(R) (in units of 
"iB/(4£p hase )) as determined numerically from Eq. (C17), vs 
R (in units of V2£ phas( ,). 

Finally, the real-space relation (C19) can be transposed 
to a local condition in wave-vector space as follows. Tak- 
ing q//z and A(q)//y (such that q • A(q) = in the 
London gauge), we obtain: 



iy(q) 



47T 



K(q)A y (q) 



where 



and 

Qyy{q) 



m B c 



m B cJ 



n B 



367r 2 e 3 



phase 



dr 



gn( r ) 1 9 / sin qr 



r 2 q dq \ qr 



(C21) 

(C22) 
. (C23) 



Following Pippard's treatment 2 ' 26 , we can then define the 
penetration depth via the relation: 



9 r°° 

A=- 



1 



(C24) 



This expression has to be compared with London's rela- 
tion at zero temperature 



Al = 



m B c z 



4irn B e B V 47rne 2 



(C25) 



(where the electric charge e B = 2e has been restored 
for convenience). To this end, we follow an analogous 
treatment within BCS theory 2 and consider the local 
(A » Cphasc) and nonlocal (A <C £ p hase) limits sep- 
arately, by analyzing the g-dependence of the kernel 
(C22). Note that in the present case, while the ratio 
Qyy(q)/Qyy(0) is a universal function of q = qV^tiphusc, 
the ratio K(q)/K(0) is not, since its asymptotic limit 
K(oo)/K(0) = 1 - (367r 2 n s ^ haS e)- 1 (with K(0) = A Z 2 ) 
is nonvanishing. This contrasts the situation in BCS the- 
ory where K(po) = 0, the difference originating from the 
fact that in Bogoliubov theory cancellation of the dia- 
magnetic term in Eq. (C22) does not occur. In the local 
(A » Cphase) limit, whereby we may approximate K(q) 
by K(0) in Eq. (C24), we obtain A ~ X L . In the non- 
local (A Cphasc) limit, on the other hand, whereby we 
may approximate K(q) by if (oo) in Eq. (C24), we obtain 
A ~ X L (K(oo)/K(0))~ 1 ^ 2 ~ \ L owing to the diluteness 
condition 46 . Within the Bogoliubov theory, there is thus 
no way to get from Eq. (C24) values for A significantly 
different from London's value A^, in contrast with BCS 
theory 2 whereby one obtains A > Al in the nonlocal 
limit. 



APPENDIX D: EXTRACTING THE PIPPARD 
KERNEL FOR SUPERCONDUCTING 
FERMIONS WITHIN THE BCS 
APPROXMATION 

In this Appendix, we set up a general procedure to ex- 
tract the Pippard kernel from the BCS bubble contribu- 
tion to the current correlation function. Contrary to the 
method introduced by BCS theory to deal with the (ex- 
treme) weak-coupling limit 2 , our procedure is applicable 
to all couplings. Although some expressions reported in 
this Appendix are similar to those of Appendix C for the 
bosonic counterpart, we shall report them here to point 
out similarities and differences between the two cases. 

In the BCS method 2 , a local (delta- function) contri- 
bution is identified from the paramagnetic contribution 
to the current induced by a (transverse) vector potential. 
This local contribution is shown to cancel completely the 
diamagnetic term, thus making "nonlocal" the relation 
between the current and the vector potential in coordi- 
nate space (over the scale of the Pippard coherence length 
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£o)- We will show in the following that this cancellation 
of the diamagnetic term by the local term originating 
from the paramagnetic contribution is, strictly speaking, 
never complete. Yet, it can be considered approximately 
complete insofar as the length scale hp 1 is much smaller 
than £o (a situation which certainly applies to the ex- 
treme weak-coupling limit). At stronger coupling (when 
kp 1 becomes comparable with £ )> on the other hand, 
this argument is bound to fail. As a consequence, the 
relation between the current and the vector potential in 
coordinate space contains a nonlocal as well as a local 
contribution. 

Quite generally, the current correlation function at zero 
frequency obtained from the BCS bubble can be cast in 
the form: 



x* c Y s (Q,n„ = o) 



i 



(2m)* J (2tt)3 (2 ^ + <W 
G BCS (q,q + Q) 



x (2qy +Qj')- 



2|q||q+Q| 



(Dl) 



where (cf., e.g., Ref. 2) 



{q,P) ~ qP E{q)E{p){E{q) + E(p)) ■ [m) 

Note that, similarly to the function G B (q,p) of Eq. (C7), 
the function G BC (q,p) is also symmetric under the in- 
terchange of p and q, vanishes for p = q, and is elsewhere 
positive. 

In coordinate space we write correspondingly 
[cf. Eq. (C6)]: 

poo poo 
-— / dqq / dppG BCS (q,p) 
V 71 "; Jo Jo 



d d 



-io^R'-RDioblR' + RI). (D3) 



2m' 

x lim 

R'^o 8W, dR' y , 

As in Appendix C, for the purpose of extracting the 
leading behavior of the function (D3) for large values 
of R = |R|, we can make the approximation (C9) in 
Eq. (D3). 

Introducing further the notation [cf. Eq. (CIO)]: 

/>oo />oo 

g BCS (R)= dqq dppG BCS \q,p) CO S (q - p)R, 
Jo Jo 

(D4) 

Eq. (D3) acquires eventually the form [cf. Eq. (Cll)]: 



_dQ 

(27T) 



3 e^ X ^(Q,O„ = 0) 



(27r) 4 m 2 i? 4 
It should be noted that also in this case 



(D5) 



dRg BCS (R) 



Jo 



dqq 2 G BCb {q,q) = (D6) 



since G BCS (q 1 q) = for all values of q. 

The integral in Eq. (D4) must be treated with care, 
since the product pG BCS (q,p) docs not vanish in the 
limit of large p, thus implying a delta- function contri- 
bution to the function g BCS (R) as for the boson case 
of Appendix C. To extract this contribution, we exploit 
the asymptotic form of the function G BCS (q,p) and set 
[cf. Eq. (C15)] 



BCS 



(q) = lim pG 

p— >oc 



BCS 



(q,p) = 8ttW 



(D7) 



where v\ = (1 - £(q)/E(q))/2 is the usual BCS factor. 
With the further notation 

G BCS (q,p) ^qG BCS (q,p)p-qf BCS ( q )-pf BCS (p), 

(D8) 



Eq. (D4) becomes [cf. Eq. (C13)]: 

g BCS (R) = g BCS (R) - 21 BCS 5(R) 
where [cf. Eq. (C17)] 



(D9) 



oo poo 



g^ b (R)= dq dpG BCb (q,p)co S ( q -p)R 
Jo Jo 

2 f°° 

+ X dqqf(q)sm(qR) (D10) 



and 



/>oo 

J BCS = -it / dqqf BCS (q) = -87r 3 mn, (Dll) 
Jo 

n being the particle density. Note that I BCS is a con- 
stant, while its bosonic counterpart of Eq. (CI 8) depend 
on the bosonic interaction. In terms of the above quan- 
tities, the relation between the induced current j and the 
vector potential A becomes eventually: 



+ 



1 

mc 
1 



jBCS 

n+ 12^)^ (r) 



(D12) 



(27r) 4 m 2 c 



^, R fY,y Ry'Aj'jr') BCS 
dr ^3 9n ( R ) 



in the place of Eq. (CI 9) for bosons. Insertion of the 
value (Dll) into the first term on the right-hand side of 
Eq. (D12) does not lead to a complete cancellation of 
the diamagnetic term, since n + T BCS /(12ir 3 m) = n/3 
in this case. In the standard BCS treatment 2 , on the 
other hand, it is claimed that the diamagnetic term is 
completely cancelled by the local (delta-function) con- 
tribution originating from the paramagnetic term. This 
apparent contradiction is resolved by the following argu- 
ment. 
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A crucial point in our analysis is that g BCS (R) is a 
regular function of R, in the sense that it does not con- 
tribute any additional delta-function term to the function 
g BCS (R), besides the one explicitly indicated in Eq. (D9). 
By a careful analysis of the one- and two-dimensional in- 
tegrals entering the definition (D10), it can, in fact, be 
shown that g„ (R) contains at most the following (in- 
tegrable) logarithmic singularity 

^ cs (i2)«-16(7r-l)m 3 A 2 In(^-) (D13) 

for Rkp <C 1. Note that the form of this singularity is 
universal, in the sense that it depends on the coupling 
value only through the value of the gap function A in the 
prefactor. 

To the singular term (D13) one should add nonsingular 
contributions that depend on the value of the coupling 
parameter. In particular, for Rkp <C 1 in weak coupling 
one obtains the nonsingular contribution Amkp. This im- 
plies that, in weak coupling, the appearance of the singu- 
lar term (D13) is restricted to extremely small values of 
R, and it cannot be detected for all practical purposes. 
For instance, when kp£o = kp/(-irmA) = 32 (so that the 
extreme weak-coupling limit is definitively reached), one 
obtains that the singular term (D13) dominates over the 
constant contribution 4mk F only for fcpi? < cxp(— 1000). 
At larger couplings, however, the singular term (D13) 
becomes non negligible for appreciable values of kpR, as 
seen in Fig. 7 of the text. [We prefer in this Appendix 
to use the Pippard coherence length £ instead of £ pa ; r 
to make easier the comparison with the BCS method to 
obtain the Pippard kernel in the extreme weak-coupling 
limit.] 



shown in Fig. 12 (full line) for kp^o = 32, whereby the 
singular behavior (D13) is not visible by the above ar- 
gument. From this plot, the two characteristic length 
scales of BCS theory (namely, kp 1 and £o such that 
kp 1 <C Co i n the weak-coupling limit) are evident. Specif- 
ically, from the plot of g BCS (R) one can isolate its "tail" 
represented by the function — (2TO7rfcF) 2 Aexp(— i?/£o) 
(or, better, by the function J(R) of BCS theory ob- 
tained by integrating the Bessel function of imaginary 
argument 2 ), shown in Fig. 12 by the dashed line. Once 
this tail has been isolated, the remainder g r l em (R) = 
9n{R) + (2mirkp) 2 Aexp(— i?/£o) is localized on the scale 
kp , which is much smaller than the scale £o of the tail 
in the weak-coupling limit. For these reasons, the re- 
mainder itself can be effectively assimilated to a delta- 
function contribution with finite strength (by setting 
g r n em {R) = -2X'S(R)), which thus adds up to the true 
delta- function contribution of Eq. (D9). Note that, to 
cancel completely the diamagnetic term in Eq. (D12), 
the factor X' should equal —Air^mn. This value is indeed 
obtained by the data of Fig. 12 within a 10~ 3 relative 
error. 

Our procedure has to be contrasted with the standard 
BCS treatment in the weak-coupling limit 2 , whereby the 
delta-function contribution that cancels completely the 
diamagnetic term in the j vs A relation is obtained by 
letting the lower limit of an integration over £ (p) and £ (q) 
going to — oo (instead of holding it at the original finite 
value —k F /(2m)). 

Although the two procedures yield correctly the same 
results in the weak-coupling limit, our procedure can 
be readily extended to the intermediate-coupling region 
where the two length scales kp 1 and £ becomes compa- 
rable, as shown in Section IV. 




5 10 15 20 25 30 35 40 45 
k F R 

FIG. 12. Plot of the function g^ cs (R) (in units of 
mk%) as determined numerically from Eq. (D4), vs kpR 
for £;f£o = 32 (full line). The corresponding func- 
tion — (2m7rfcF) 2 Aexp(— R/£o) is also shown for comparison 
(dashed line). 

A plot of the function (R) (in units of mk F ) vs 
kpR obtained by numerical evaluation of Eq. (D10) is 
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